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Abstract 

We conjecture that the structure of Bernoulli numbers can be explicitly given in 
the closed form 

s n =(-i)t- 1 n M" 1 n ipcxwj-ftjip 1 n^ 1 

p— X\n (p,0e*i rr p-l\n 

n = l (modp-1) 

where the X(p,l) are zeros of certain p-adic zeta functions and ^ff 1 is the set of irregular 
pairs. The more complicated but improbable case where the conjecture does not hold 
is also handled; we obtain an unconditional structural formula for Bernoulli numbers. 
Finally, applications are given which are related to classical results. 
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1 Introduction 

The classical Bernoulli numbers B n are defined by the power series 

oo n 

— 7=y, B n— 7, \Z <27T, 



n=0 



where all numbers B n are zero with odd index n > 1. The even- indexed rational 
numbers B n alternate in sign. First values are given by Bq = 1, B\ = — |, B2 = |, 
B4 = — ^j. Although the first numbers are small with \B n \ < 1 for n = 2,4, . . . , 12, 
these numbers grow very rapidly with \B n \ — > 00 for even n — > 00. 

For now, let n be an even positive integer. An elementary property of Bernoulli numbers 
is the following discovered independently by T. Clausen [Cla40] and K. G. C. von Staudt 
[vS40] in 1840. The structure of the denominator of B n is given by 

B n + ^ -eZ and denom(B n ) = [ p. (1.1) 

p— l\n ^ p— \\n 
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As a result, often attributed to J. C. Adams, see [IR90, Prop. 15.2.4, p. 238], B n /n is 
a p-integer for all primes p with p — 1 \ n. Therefore, one has a trivial divisor of B n 
which is cancelled in B n /n 

Y[ p ord P n | B n . (1.2) 

p— lfn 

This is now known as Adams' theorem, although Adams [Ada78] only predicted without 
proof that p \ n implies p \ B n for primes p — 1 \ n regarding the table of the first 62 
Bernoulli numbers i?2m he had calculated in 1878. On the other side, the property 
that B n /n is a p-integer for p — 1 \ n is necessarily needed to formulate the so-called 
Kummer congruences given by E. E. Kummer [Kum51] earlier in 1851. 

These congruences and its generalizations are important properties of Bernoulli num- 
bers which lead to a p-adic view giving interesting information about B n /n. Let ip be 
Euler's totient function. The Kummer congruences state for n,m,p,r 6 N, n,m even, 
p prime and p — 1 \ n 

(l_p«-i)^L = (l_ p n^iA r modp r^ ( 13 n 
n m 

with n = m (mod ip{p r )). 

In 1850 Kummer [Kum50] introduced the classification of regular and irregular primes 
to characterize solutions of the famous Fermat's last theorem (FLT). An odd prime p 
is called regular if p does not divide the class number of the cyclotomic field Q(/i p ) 
with /i p as the set of p-th roots of unity, otherwise irregular. Kummer proved that if 
p is regular then FLT has no solution for the exponent p. He also gave an equivalent 
definition concerning Bernoulli numbers: An odd prime p is called regular if p does not 
divide any Bernoulli number B n for n = 2,4, . . . ,p — 3, otherwise irregular. The index 
of irregularity i(p) counts these indices for which p \ B n happens. In this case the pair 
(p, n) is called an irregular pair. First irregular primes are 37, 59, 67, 101. 

Regarding Bernoulli numbers, it will be very useful to combine properties of B n as well 
of B n /n, the so-called divided Bernoulli number. An easy consequence of the Kummer 
congruences provides that the numerator of B n /n consists only of irregular primes and 
that infinitely many irregular primes exist. For the latter see a short proof of Carlitz 
[Car54], see also [IR90, Theorem 6, p. 241]. Unfortunately, the more difficult question 
is still open whether infinitely many regular primes exist. However, calculations of 
Buhler, Crandall, Ernvall, Metsankyla, and Shokrollahi [BCE + 01] show that about 
60% of all primes less 12 million are regular which agree with an expected distribution 
proposed by Siegel [Sie64]. 

All these basic results of Bernoulli numbers can be found in the book of Ireland and 
Rosen [IR90, Chapter 15]. Throughout this paper all indices concerning Bernoulli 
numbers will be even and p an odd prime. Note that in older references the enumerating 
of Bernoulli numbers can differ by a factor 2. Let p r \\ n denote the highest power of 
p dividing n in order that r = ord p ra. 
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2 Preliminaries 



Here we will recall necessary facts about irregular prime powers of Bernoulli numbers 
and p-adic zeta functions. The definition of irregular pairs can be extended to irregular 
prime powers which was already given in [Kel04], first introduced by the author [Kel02, 
Section 2.5]. 

Definition 2.1 A pair (p,l) is called an irregular pair of order n if p n \ Bi/l with 
2 < I < f(p n ) and even I. Let 

K r :={(p,l):p n \B l /l, 2<l<v(p n ), 2 | 1} 

be the set of irregular pairs of order n. For a prime p the index of irregular pairs of 
order n is defined by 

in(p) := #{(p,l): (P,0e*™}. 
Let (p, I) S ^JJ r be an irregular pair of order n. Let 

n 

( P , Si, 8 2l ■ ■ ■ , S n ) € &™ , l = J2 Sy ^{p"' 1 ) 

v=l 

be the p-adic notation of (p, I) with < s v < p for v = 1, . . . , n and 2 | si, 2 < si < p— 3. 
The corresponding set will be denoted as ^" T - The pairs (p,l) and (p, si, S2, ■ ■ ■ , s n ) 
will be called associated. Let (p, I) € ^> l " be an irregular pair of order n. Then define 

with < A(p_j) < p. In the case A( p ^ = we will denote A^n as singular. 

Remark 2.2 Note that this definition includes the usual definition of irregular pairs 
for n = 1 with i{jp) = i\(p). By Kummer congruences (1.3) the interval [2,ip(p n ) — 2] 
is given for irregular pairs of order n if they exist. Moreover, we have the property 
that if (p,l) £ then p n \ Bi +kip ^/(l + kip{p n )) is valid for all k G N . For 
simplification (p, s\,S2, ■ ■ ■ , s n ) is also called an irregular pair with (si, S2, ■ ■ ■ , s n ) as 
the second parameter in a p-adic manner. It is easy to see that if (p,l n+ i) E ^n+i 
exists then (p,l n ) 6 exists with (p, l n ) = (p, l n +i (mod </?(p n ))), too. 

The following proposition, see [Kel04, Prop. 5.3], gives an unconditional representation 
of Bernoulli numbers by sets \I/JJ r . This is seen by (1.1), (1.2), and counting irregular 
prime powers. 

Proposition 2.3 Let n be an even positive integer, then 

B n = (-l)^ 1 ft p^»)-^» / 11 p 
p— X\n p— l\n 

with 

oo 

r(p,n) := ^#(^ rr n{(p,n (mod <p(p v )))}) ■ 

v=l 
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The divided Bernoulli numbers B n /n are directly related to the Riemann zeta function 
on the negative x-axis 

C(l-n) = -, nGN,n>2, (2.1) 

n 

where the Riemann zeta function is usually defined by the sum or the Euler product 

oo 

C(s) = J2^ s = ll^-P~T\ seC,Res>l. (2.2) 

v=\ p 

Let Z p be the ring of p-adic integers and Q p be the field of p-adic numbers. Let | \ p be 
the ultrametric absolute value on Z p which is usually defined by \x\ p = p~ ord P x . The 
property of Kummer congruences allows the construction of a unique continuous p-adic 
zeta function which was introduced by T. Kubota and H. W. Leopoldt [KL64] in 1964, 
see also Koblitz [Kob96, Chapter II]. Modifying (2.1) with an Euler factor (1 — p 11 ^ 1 ) 
implies the following definition. 

Definition 2.4 Let p be a prime with p > 5. Let 

Cp(l-n) := (l-p^KCl-n) = (1-p"" 1 ) (-^) ■ 

For a fixed s± G {2, 4, . . . , p — 3}, define the p-adic zeta function by 

C P , S1 : Z p -> Z p , Cp, si (s) := lim Cp(l - (si + (p - 1)^)) 

for p-adic integer s by taking any sequence (t u ) u >x of nonnegative integers which p- 
adically converges to s. 

By construction the p-adic zeta function £ p S1 (s) interpolates the zeta function Q p (l — n) 
at nonnegative integer values s by 

( P , si (s) = C P (1 - n) 

with n = si (mod p — 1) and n = s\ + (p — l)s. Because of Kummer congruences 

Cp(l - n) = Cp(l - n') (mod p r ) 

for n = n' (mod <f(p r )) with n = re' = si (mod p — 1), the p-adic zeta function Cp,si ( s ) 
is a unique continuous function on Z p by means of interpolating property. 
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3 The A-Conjecture 



Let (p,l) € tyf 1 be an irregular pair, then we are interested in the behavior of 
Essentially, there are two cases to consider: A( p n ^ and Ar p n = 0. Now, calculations 
in [BCE+01] for irregular primes p < 12 000 000 show that A^n ^ is always valid. 
No singular A^) has been found yet. However, the improbable case of a singular A^n 
which implies a strange behavior without regularity is described in the next section. 
The following theorem gives the main result of irregular pairs in the nonsingular case, 
see [Kel04, Theorem 3.1]. 

Theorem 3.1 Let (p,h) £ ^i rr be an irregular pair with A( p ^) ^ 0. Then for each 
n > 1 there exists exactly one irregular pair of order n corresponding to (p,li). There- 
fore, a unique sequence (l n )n>i resp. (s n ) n >i exists with 

(P, In) € K T resp. (p, si, . . . , s n ) G §£ r 

and 

h < h < ^3 < • • • ■ lim l n = oo . 



A (pM) = A (P,h) = A (p,/ 3 ) 



Moreover, one has 

If Afpfa s 7^ for all i{p) irregular pairs (p, l\ jV ) 6 ^i T , then 

i(p) = i 2 (p) = h{p) = ■■■ ■ 

Definition 3.2 Let (p, I) £ fy 1 ™ be an irregular pair with ^( Pj i) ^ 0. Then Theorem 
3.1 provides a unique sequence (s u ) u >i with I = s\. Define a characteristic p-adic 
integer 

which contains all information of irregular pairs of higher order corresponding to (p, I). 

The following theorem, a result of [Kel04, Theorem 4.6/4.10], shows the behavior of 
the p-adic zeta function. The Kummer congruences (1.3) are valid by the implication 

n = m (mod <p(p r )) =^ (1 - p* 1 ' 1 )^ = (1 - p™' 1 )^ (mod/), 

n m 

but the converse does not hold in general. The first nontrivial counterexample is given 
by p = 13 and 5i 6 /16 - B 4 /A = -7 • 13 2 /2720. Note that also Si 4 /14 - B 2 /2 = 
happens which is the only exception that divided Bernoulli numbers are equal. 

Theorem 3.3 Let (p, I) € \I ,1 1 rr be an irregular pair with A( p ^ ^ 0. Let s, t 6 Z p . Then 
the p-adic zeta function has a unique zero with Cp,i{x(p,i)) = 0. ^4 strong version 
of the Kummer congruences holds 

\C P ,i(s) ~ Cp,i(t)\ p = \p(s-t)\ p . 

Moreover, one has 

( P ,i(s) - ( P ,i(t) 



p(s-t) 



-A(p,i) (mod p7L v ) , s^t. 
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As a consequence, we easily obtain |Cp,«( s )lp = \p(X(p,l) ~ s )\p f° r s ^ TL V under the 
assumption above. Thus, the power of the irregular prime p can be described by 
measuring the p-adic distance to the zero X(v,l) °f the p-adic zeta function ( p j. Since 
\pn\ p |C(1 - n)\ p = 1 for n = (mod p - 1) and |C(1 - ri)\ p = \( P ,i{s)\ p for s = 
n = / (mod p — 1), we obtain a structural formula of £(1 — n) resp. B n /n, see [Kel04, 
Theorem 4.9]. Note that we combine the numerator and denominator of B n /n in this 
formula. Vaguely speaking, the numerator can be described by zeros of £p,z and the 
denominator by poles of £ P ,o lying at 0, where ("p, o : <^p — ► Qp extends Definition 2.4 
using arguments given in [Kob96, Chapter II, p. 46]. 

Theorem 3.4 Let P be the set of primes. Then define = *i rr U (P x {0}) and 
X(p,o) = f or a M P G P- Define p(l) = 1 — 2sign(/) = ±1 for I > 0. Let n be an even 
positive integer, then under the assumption that no singular n exists, one has 

C(l-n) = (-l)? n 

n = l (modp— 1) 

One may ask whether the structure of the Riemann zeta function at odd negative 
integers is given by this simple form. Now, all these facts substantiated by calculations 
lead to the following conjecture, already mentioned in [Kel02]. 

Conjecture 3.5 (A-Conjecture) Let p be an irregular prime. Then the following 
properties hold: 

(1) A(p ^) is not singular for all irregular pairs (p,l u ) £ ^f 1 - 

(2) i{p) = i 2 (p) = h{p) = .... 

(3) The p-adic zeta function Cp,i„ has a unique zero X(p./„)- 

(4) A strong form of the Kummer congruences holds for Cp,i v - 

Finally, assuming the A-Conjecture, we also obtain a structural formula of Bernoulli 
numbers which gives a decomposition in three products. The first and last product are 
trivially given, the complicated product in the middle consists only of irregular primes. 

Theorem 3.6 Let n be an even positive integer, then under the assumption of the 
A-Conjecture 

B n =(-i)^- 1 n hj 1 n iPixM-^ 1 Up 1 

p— l\n (p,Z)e*Y" r p—l\n 

n = l (modp— 1) 

Proof. The product in the middle is given by Theorem 3.4. The first resp. last 
product is a consequence of (1.2) resp. (1.1). □ 

The main task remains to determine the zero of a p-adic zeta function associated with an 
irregular pair (p, /). An irregular pair (p, l n ) G of order n yields an approximation of 
the zero X(p, i)- Fortunately, these irregular pairs of higher order can be computed with 



|X(p,0 



l lp 



P 
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little effort by calculating a small number of divided Bernoulli numbers with relatively 
small indices. For algorithms and calculated pairs see [Kel04, Section 5, Table A. 3]. For 
example, we have 

(157, 62, 40, 145, 67, 29, 69, 0, 87, 89, 21) € 

which also shows a rare occurrence of a zero in the p-adic sequence, here sj = 0. This 
is the only zero which occurs in the p-adic sequence of irregular pairs of order 10 for 
p < 1000, see [Kel04, Table A.3]. This means 

(157, 6 557 686 520 486) e &™ n # 7 rr 

respectively, since the index is prime to 157 and sg = 87, 

157 7 || Bq 557 6 86 520 486 • 

Now, we will never be able to compute this giant Bernoulli number nor we can use 
Kummer congruences to verify this result directly! 



4 The singular case 

The case of a singular A^n is more complicated than the nonsingular case where 
we have a certain regularity. However, no such singular A^n has been found yet. 
Therefore, supported by massive computations, we can regard such an event as a rare 
exception. The following theorem, see [Kel04, Theorem 3.2], shows the strange behavior 
of irregular pairs of higher order with a singular A( p n. 

Theorem 4.1 Let (p, l n ) € be an irregular pair of order n with A(p t i n ) = 0. Then 
there exist two cases: 

(1) (p,l n ) £ ^^4-1 : There are no irregular pairs of order n + 1 and higher. 

(2) (p,l n ) G *n+i : There exist p irregular pairs (p,l n+1 j) = (p,l n + j<p(p n )) G ^"+1 
of order n + l with A {P)/n+1 ^ = for j = 0, . . . ,p - 1. 

The following diagram demonstrates the behavior of the singular case. This situation 
can be described by a rooted p-ary tree of irregular pairs of higher order. 

w 
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Here a vertical line indicates that (p, l n ) G ty 1 ™ H happens. We then have p 

irregular pairs of order n + 1 which are represented by branches. In this case, the 
corresponding Bernoulli number Bi n /l n decides whether there exist further branches 
or they stop. Instead of n the order of the p-power must be at least n + 1. This 
also means that an associated irregular pair (p, si, . . . , s n +i) 6 must have a zero 

s n +i = in its p-adic notation each time. Now, it is worth saying that no irregular 
pair (p,l) has been found with p 2 \ B t resp. (p,Z) € *f r D *£ r for p < 12 000 000, 
see [BCE + 01], while an example of an element of \I/ 6 rr n ^>" 1 is shown in the previous 
section. 

Definition 4.2 Let (p, I) G \E ,1 1 rr be an irregular pair with a singular A( p n. Then define 
a rooted p-ary tree of irregular pairs of higher order like in the diagram above given by 
Theorem 4.1. Each node contains one irregular pair of higher order. Note that these 
pairs are not necessarily distinct. We denote this tree as T? « related to the root node 
(p, I). The tree T^ p ^ has the property that each node of height r lies in Vf^+i- ^ * ree 
T9 « = {(p, /)} is called a trivial tree having height 0. A tree of height one is given by 
the root node (p, I) and its p child nodes (p, I + jip(p)) with j = 0, . . . ,p — 1. A tree 
with height > 2 always contains the latter one. 

In the nonsingular case, we have a zero of the p-adic zeta function. In contrast to, 
the singular case does not guarantee that irregular pairs of higher order exist at all. 
However, an exception does not destroy Theorem 3.6 but complicates the formula, 
because we then have to consider the tree T? « of irregular pairs of higher order. 
Thus, we obtain an unconditional formula by combining both cases. 

Theorem 4.3 Let n be an even positive integer, then 

B n = (-i)t- 1 n M- 1 n [p( X(P)0 -^)|-i 

p-l\n (p,i)e*j r , A d»,!)^ 

n=l (modp — 1) 

p l+h (p,n) J-j- p -X 
( P ,Oe*i",A (piir o p-l\n 

n=l (modp— 1) 

with the height of (p, n) defined by 

h (p,n) = max<jheight((p,0) : (p, € T^ l} n {(p,n (mod p(p")))}„>i} . 
Moreover, ho(p,n) = i/ie £ree is trivial. Otherwise ho(p,n) > 1. 

Proof. The case A( p ^ 7^ is already handled by Theorem 3.6. Now, assume A(p n = 
with a given tree T^iy We have to determine the max. height of a node (p,l u j) G 
TjP z) ^ wn i c h is equal to (p, n (mod (^(p^))), as a consequence of Remark 2.2. The 
root node (p, I) has height 0, so the exponent equals 1 + ho(p,n). If the tree TtLn is 
trivial, then ho(p,n) = is constant. On the other side, a tree having height > 1 
contains p irregular pairs of order two. Then (p, n (mod (p(p 2 ))) G TP ^ is always valid 
which finally yields ho (p, n) > 1. □ 
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Corollary 4.4 Let n be an even positive integer, then 

ca-n) =(-i)§ n ipcx^o-^Ip 1 n ^ i+Mp,n) n ^ 

(p,i)e*i rr .A (pii) #o ( P ,oe*i rr ,A (p0 =o p-l|n 

n = i (mod p — 1) n = i (mod p— 1) 

wii/i ho(p,n) as defined above. 

5 Applications 

Regarding Theorem 4.3 and the definitions of ho and X(p,0' we can state an extended 
version of Adams' theorem given by (1.2). 

Theorem 5.1 Let n be an even positive integer. Let p be a prime with p r \ \ n, r > 1, 
and p — 1 \ n. Let I = n (mod p — 1) with < Z < p — 1. TZien p r+<5 || B n with the 
following cases: 

(1) If p is regular, then 5 = 0. 

(2) If p is irregular with (p,l) ^ Vf^, then 5 = 0. 

(3) Ifp is irregular with (p, I) G *j r ; A (Pi0 ^ 0, i/ien 5=1 + ord p (x(p,j) - ^Er). 

(^j //p is irregular with (p,l) G \I /1 1 rr ; = 0, then 5=1 + ho(p,n). 

Additionally, in case (3) resp. (4), if (p,l,l) £ then 5 = 1 is bounded, otherwise 
5>2. 

Proof. We have to consider the formula of Theorem 4.3, then the first product yields 
p r | B n . Only the second resp. third product can give additional p-factors. Therefore, 
case (1) and (2) are given by definition. Now, we can assume (p, I) G ^f 1 . A nonsingular 
A(p,/) provides 5 = ord p \p (x( P ,i) ~ p*ff)lp 1 = 1 + °rd p {x(p,l) ~ f=r) in case (3). On the 
other side, a singular A^n provides 5 = 1 + ho(p, n) in case (4). The additional cases 
are shown as follows. 

Case (3): By assumption, n = p r n' with some integer n'. We have to evaluate 

d = ordp (x(p,l) ~ %Ej) = ord P (PX(p,0 - X( P ,l) + 1 ~ tfn) . 

Since r > 1, we p-adically obtain X(p,l) = I + s^p + . . . -4=>- d > 1 which is equal to 
(p,l,l) G % r . Conversely, (p,l,l) £ % r yields d = 0. Case (4): (p,l,l) £ % r 
the tree T9 ^ is trivial. Then ho(p, n) = is constant. Conversely, (p, I, I) G VPJf 1, yields 
ho(p, n) > 1 as a result of Theorem 4.3. □ 

So far, no (p, Z, Z) G Vl/Jp has been found yet. We can even raise the value 5 in the 
following way. 

Corollary 5.2 Assume that (p,l, . . . ,Z) G ^JT 1 exists with some r > 1. ZeZ n = Zp r . 
Then, we have p r \ \ n and p 2r \ B n . 
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Proof. By Definition 2.1, we have 



i r = j2 1 ^(p^ 1 ) = ¥~ l with ( P , i r ) g . 

Then p r \ B lr+klfi ( p r}/(l r + k (f(p r )) is valid for all k > 0. Choose n = l r + l (p(p r ) = lp r . 
Thus, p r | B n /n and finally p 2r | S n . Note that we cannot predict that p 2r \\ B n in 
general. □ 

Johnson [Joh74, Theorem, p. 655] calculated irregular pairs up to p < 8000. Corre- 
spondingly, he also calculated the now called irregular pairs (p, s±,S2) € ^> 2 " of order 
two in that range, proving that (p, I, I) ^ ty 1 ™ f° r P < 8000. See also [Kel04, Table A. 3] 
for calculations of irregular pairs of order 10 for p < 1000. In a similar manner, the 
nonexistence of irregular pairs (p, I, I — 1) of order two plays an important role in Iwa- 
sawa theory, see Washington [Was97] for Iwasawa theory and [Kel04, Section 6] for this 
special result. In context of cyclotomic invariants, calculations of [BCE + 01] ensure that 
no (p, 1,1 — 1) £ exists for p < 12 000 000. One may conjecture that no such special 
irregular pairs (p, /, I) resp. (p, 1,1 — 1) of order two exist. But there is still a long way 
to prove such results, even to understand properly which role the zeros X(p,i) play- 

By Definition 2.1, we have the relation 

(p,l,l) i W p 2 \B lp /{lp) <=>> p*\B lp . 

Yamaguchi [Yam76] also verified by calculation that p 3 \ Bi p for all irregular pairs 
(p, /) with p < 5500, noting that this was conjectured earlier by Morishima in general. 
Interestingly, the condition p 3 { B[ v is related to the second case of FLT, see [Was97, 
Theorem 9.4, p. 174]. See also [Was97, Corollary 8.23, p. 162] for a different context. 
Under the assumption of the conjecture of Kummer-Vandiver and that no (p, 1,1) & 
exists, the second case of FLT has no solution for the exponent p. For details we refer 
to the references cited above. 

Remark 5.3 Thangadurai [Tha04] also investigates Adams' theorem claiming a con- 
jecture that 5 € {0,1} where 5 is given as in Theorem 5.1. The following theorem 
[Tha04, Theorem 2.5, p. 172] is formulated, here translated into our terminology 

5 G {0, 1} <^> p 3 \ B ip p 2 \ B\ for all I = 2, 4, . . . ,p - 3 . 

The first equivalence agrees with our results, but the second equivalence is false. The- 
orem 2.3 and Theorem 2.4 of [Tha04, pp. 171-172] are incorrect in general. Theorem 
2.3 of [Tha04] is an incorrect citation of the results of Johnson [Joh74, Theorem, p. 655] 
which are only valid for irregular primes p < 8000 and were verified by calculations. 
In our words, Theorem 2.3 of [Tha04] would imply A( p ^ / => (p, I, I) ^I' 1 2 rr for 
all irregular pairs! But A( p ^ ^ only implies the existence of a (p,l,k) S with 
a unique k in the range < k < p. Theorem 2.4 of [Tha04] states the following: Let 
(p,l) € ^i rr . Ifp || B l+k (p_^/(l + k(p-l)) for some k > 1, thenp || B L . Now, this theo- 
rem would imply that p 2 \ B[ in the nonsingular case! Counterexample: A( p n ^ ==r- 
there exists exactly one (p,l,k f ) G ^f 1 "- The case k' = can occur. Then we have 
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p 2 I Bi, but p || Bi + Lp_i^/{1 + (p — 1)). As a consequence, [Tha04, Theorem 1, p. 170] 
is only valid for irregular primes p < 8000 by results of Johnson. Finally, [Tha04, The- 
orem 4.2, p. 176] is also incorrect, because (p, 1,1 — 1) ^ VPip does not imply p 2 \ Bi in 
the nonsingular case A( p n 7^ 0, while in the singular case this implication always holds 
as seen by a trivial tree TfLiy I n conclusion, all depending results of these theorems 
cited above are incorrect, especially the claimed equivalences to other hypotheses of 
Bernoulli numbers in [Tha04]. 

Certainly, the converse of Adams' theorem does not hold, but one can state a some- 
what different result which deals with the common prime factors of numerators and 
denominators of adjoining Bernoulli numbers, see [Kel02, Satz 2.3.4, p. 35]. 

Theorem 5.4 Let S = {2, 4, 6, 8, 10, 14} be the set of all even indices m where the 
numerator of B m /m equals 1. Write B n = A n /V n with (A n ,V n ) = 1. Let k,n be even 
positive integers with k G S and n — k > 2. Then 

D = (A n ,V n -k) provides D \ n . 

Moreover, if D > 1 then D = p\ ■ ■ ■ p r with r > 1 and p v \ V\-, p v \ B n /n for all 
v = l,...,r. 

Proof. Assume D > 1. We then have D = p\ ■ ■ ■ p r with r > 1 since V n -t is squarefree 
by (1.1). Let v £ {1, . . . ,r}. We have the following properties: p u \ A n and p v — 1 { n, 
additionally, p u \ V n -k and p v — 1 | n — k with p u < n. Hence, we obtain p u \ Vfc. 
Assume to the contrary that p u — l\k. By p u — 1 | n — k we get p u — 1 | k + n — k = n. 
Contradiction. 

Assume p v \ n or p u \ B n /n, then we obtain by Kummer congruences (1.3) 

B n Bk , . 
0= — =— (modp^J, 
n k 

since n — k = (mod p u — 1). By properties of the set S 

-^0 (modp u ) (5.1) 

yields a contradiction. We obtain p u | n and p v \ B n /n. Finally D \ n is valid. □ 

Now, the set S cannot be enlarged, because (5.1) does not hold in general for numera- 
tors having prime factors. For example, let p = 691 and n = 12 + {p — 1), then we have 
p I B12/I2 and D = (A n , V n -i2) = pc\ n with some c > 1. Actually, c = 1 with the 
help of Mathematica. On the other hand, one trivially obtains for k G S, p prime 
with p — 1 \ k, n = kp infinitely many examples of D > 1. In the following theorem, 
Theorem 5.4 plays a crucial role. Define for positive integers n and m the summation 
formula of consecutive integer powers by 

m— 1 

S n (m) = J2» n - 

Many congruences concerning function S n are naturally related to Bernoulli numbers. 
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Theorem 5.5 Let n, m be positive integers with even n. Then 



m 



r+l 



S n (m) 



m 



B r 



for r 



1, 2. 



PROOF. Write B n = A n /V n with (A n ,V n ) = 1. Assume m > 1 and n > 10 with even 
n, otherwise we have \A n \ = 1 for n = 2,4,6,8. It is well known, see [IR90, p. 234], 
that for even n > 10 

S n (m) =B n m+ (^J B n „ 2 — + £ (^J ^n-fc-^y • (5.2) 

We have to examine carefully the sum given in (5.2). By the result of Clausen-von 
Staudt (1.1) the denominator of all nonzero Bernoulli numbers is squarefree including 
So an d B\. Consider each prime factor p s \ \ m. Then, we have 

{ fn\ m k+1 \ [ s > k > 2 >P> 2 

ordpM k ) B n-hj^J >s{k + l)-l-ord p {k + l)> I 2s,k>2,p>5 (5.3) 

I 3s , k ^ 4, p ^ 5 

for such k < n where .B n _fc ^ 0. Critical cases are to consider for p = 2, 3, 5 and s = 1. 
Now, we are ready to evaluate (5.2) (mod m') for certain t. 

Case r = 1: Assume (m, V n ) > 1, then we obtain by (5.3) (case k > 2, p > 2) 

S n (m) = B n m = — r m ^ (mod m) . 

Therefore, (m,V n ) = 1, 2 | m, 3 { m, and p > 5 must hold. Hence, by (5.3) (case 
k > 2, p > 5), we can write S n (m) = B n m (mod m 2 ). Consequently, 

= S n (m) = S n m (mod m 2 ) 

provides m 2 | S n (m) <^=^> ?n | B n . 

Case r = 2: We have m \ B n and (m, 6) = 1, because either m 2 | B n or m 3 | S n (m) 
is assumed. The latter case implies m 2 \ S n (m) and therefore with case r = 1 also 
m | B n . By (5.3) (case A; > 4, p > 5) we obtain 

n\ B n ^ 2 _ 3 _ D n(n - l)vl n _ 2 3 



S n {m) = B n m+\ \^-^m 6 = B n m+ v y m 3 (mod m 3 ) . (5.4) 
\2/ 3 6V/ n _ 2 

Our goal is to show that the second term vanishes, but the denominator V n ~2 could 
possibly remove prime factors from m. Now, Theorem 5.4 asserts that (A n ,V n -2) \ n. 
We also have (m, V n -2) \ n since m \ B n . This means that the factor n adds those 
primes which possibly removes from m. Therefore, the second term of (5.4) 

vanishes (mod m 3 ). The rest follows again by S n (m) = B n m = (mod m 3 ). □ 

One can improve the value r for certain m \ n, since (^) appears in the sum (5.2), but 
not in general. Let p = 37 and / = 37580. We then have (p, I) € VPlp and p 3 | Bi, but 
p 4 \ Si(p) which was checked with Mathematica. 
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Example: 

(1) We have B 42 = 1520097643918070802691/1806. Since the numerator A 42 is a big 
prime, we obtain for m > 1 

m 2 | S i2 (m) ^> m = 1520097643918070802691 . 

(2) We have A 5Q = 5 2 • 417202699 • 47464429777438199 and V 48 = 2 • 3 • 5 • 7 • 13 • 17. 
Hence, for m > 1 

m 3 | S l 5o(m) m = 5 . 

Bernd C. Kellner 

address: Reitstallstr. 7, 37073 Gottingen, Germany 
email: bk@bernoulli.org 
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